Microscopic theory of quantum dot interactions with quantum hght: local field effect 
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A ttieory of botli linear and nonlinear electromagnetic response of a single QD exposed to quantum 
light, accounting the depolarization induced local-field has been developed. Based on the micro- 
scopic Hamiltonian accounting for the electron-hole exchange interaction, an effective two-body 
Hamiltonian has been derived and expressed in terms of the incident electric field, with a separate 
term describing the QD depolarization. The quantum equations of motion have been formulated 
and solved with the Hamiltonian for various types of the QD excitation, such as Fock qubit, coherent 
fields, vacuum state of electromagnetic field and light with arbitrary photonic state distribution. For 
a QD exposed to coherent light, we predict the appearance of two oscillatory regimes in the Rabi 
effect separated by the bifurcation. In the first regime, the standard collapse-revivals phenomenon 
do not reveal itself and the QD population inversion is found to be negative, while in the second 
one, the collapse-revivals picture is found to be strongly distorted as compared with that predicted 
by the standard Jaynes-Cummings model. For the case of QD interaction with arbitrary quantum 
light state in the linear regime, it has been shown that the local field induce a fine structure of the 
absorbtion spectrum. Instead of a single line with frequency corresponding to which the exciton 
transition frequency, a duplet is appeared with one component shifted by the amount of the local 
field coupling parameter. It has been demonstrated the strong light -mater coupling regime arises 
in the weak-field limit. A physical interpretation of the predicted effects has been proposed. 



PACS numbers: 42.50.Ct,73.21.-b,78.67.Hc 



I. INTRODUCTION 

The strong coupling between condensed matter and 
quantum light is a core issue of present day quan- 
tum optics. Realized by exposing the matter with in- 
tense quantum field, it can manifest itself in differ- 
ent quantum systems such as single atoms and ultra- 
cold atomic beamsi, semiconductor hctcrostructures2^3, 
Bose-Einstein condensates^, etc. Albeit these systems 
are of different physical nature, their interaction with 
quantum light is governed by common rules. In the 
strong coupling regime, these systems enable the gen- 
eration of different states of quantum light — single 
photons^, Fock states^, Fock qubits, quantum states with 
arbitrary photon number distribution^. That constitutes 
a basis for the quantum information processing^i^ and 
quantum metrology^*^. In practice, the strong light- 
matter coupling regime can be realized in two ways: by 
combining the matter with a high-Q microcavity or by 
exposing the matter to a ultrashort intense pump pulse. 

To describe the strong coupling between an arbi- 
trary two-level system and quantum light, the Jaynes- 
Cummings (JC) model is conventionally usedii. One 
of the most fundamental phenomenon predicted within 
the JC model is the oscillation of the population be- 
tween levels with the Rabi frequency (Rabi oscillations). 
However, the standard JC model does not account 
for a number of physical factors, which, under certain 
conditions, may significantly influence the Rabi effect. 
The time-domain modulation of the field-matter cou- 
pling constant-^J^ and interplay between classical driv- 



ing field and quantized cavity field^ can serve as exam- 
ples. More advanced JC models involve additional inter- 
action mechanisms and effects, such as dipole-dipole (d- 
d) interactior^i^iii, exction-phonon coupulingi^ii^, and 
self-induced transparency^. The d-d interaction be- 
tween two quantum oscillators leads to radiative coupling 
of them and, as a result, to exchange by the excited state. 
That is, Rabi oscillations between these two oscillators 
occur; see Ref. [H for a theory and Ref. [l^ for the exper- 
imental observation in a double quantum dot (QD) sys- 
tem. As a whole, the observation and intensive studying 
of excitonic Rabi oscillation a^^i^°i^-^i^^i^^i^'^i^^ motivates 
the extension of the JC-model to incorporate specific in- 
teractions inherent to confined exciton in a host. 

In the given paper we present a microscopic theory of 
the interaction of an isolated QD with quantum light for 
both weak and strong coupling regimes. We incorporate 
the local field correction into the JC model as an addi- 
tional physical mechanism influencing the Rabi effect in 
a QD exposed to quantum light. In particular, the Rabi 
oscillations are shown to exist even in the limit of a weak 
incident field. 

In the weak coupling regime, the local-field effects in 
optical properties of QDs have been theoretically inves- 
tigated in Refs. |2ll23!28'^29',30'3T'35'33"3i for classical 
exposing light and in Ref. 35 for quantum light. In the 
latter case it has been shown that for a QD interacting 
with Fock qubits the local fields induce a fine structure of 
the absorption (emission) spectrum: instead of a single 
line with the frequency corresponding to the exciton tran- 
sition, a doublet appears with one component shifted to 
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the blue (red). The intensities of components are com- 
pletely determined by the quantum light statistics. In 
the limiting cases of classical light and Fock states the 
doublet is reduced to a singlet shifted in the former case 
and unshifted in the latter one. 

The role of local fields in the excitonic Rabi oscilla- 
tions in an isolated QD driven by classical excitation 
was investigated in Ref. i36j. Two different oscillatory 
regimes separated by the bifurcation have been predicted 
to exist. The Rabi oscillations were predicted to be non- 
isochronous and arising in the weak excitation regime. 
Both peculiarities have been experimentally observed by 
Mitsumori et al. in Ref. [2^ where the Rabi oscillations of 
excitons localized to quantum islands in a single quantum 
well were investigated. 

There exist several different physical interpretations of 
local field in QDs and, correspondingly, different ways of 
its theoretical description. The first model (scheme A in 
the terminology of Ref. Is^) exploits the standard elec- 
trodynamical picture: by virtue of external field screen- 
ing by charges induced on the QD surface (the qua- 
sistatic Coulomb electron-hole interactions) a depolar- 
ization field is formed differentiating the local (acting) 
field in the QD and external incident field. In this model 
the total electromagnetic field is not pure transverse. Al- 
ternatively, only transverse component is attributed to 
the electromagnetic field, while the longitudinal compo- 
nent is accounted for through the exchange electron-hole 
interactions (scheme B accordingly Ref. Isol ). Both ap- 
proaches are physically equivalent and lead to identical 
results. 

In the present paper we build the analysis on the gen- 
eral microscopic quantum electrodynamical (QED) ap- 
proach where the local field correction originates from the 
exchange by virtual vacuum photons between electrons 
and holes forming the exciton and thus is a manifestation 
of the dipole-dipole (d-d) interaction between electrons 
and holes (the dynamical Coulomb interaction) The 
approach allows us to overcome a number of principal 
difficulties related to the field quantization in QDs^. 
In the analysis, approximate solution of the many-body 
problem is built on the Hartree-Fock-Bogoliubov self- 
consistent field concept The self-consistent technique 
leads to a separate term in the effective Hamiltonian re- 
sponsible for the interaction of operators and average val- 
ues of physical quantities. Due to this term the quantum 
mechanical equations of motion become nonlinear and 
require numerical integration. 

The paper is arranged as follows. In Sect. |TT] we de- 
velop theoretical model describing the QD-quantum light 
interaction. We formulate a model Hamiltonian with the 
separate term accounting for the local field correction and 
corresponding equations of motion. In Sect. IIIII we ana- 
lyze the manifestation of local fields in the motion of the 
QD exciton in the absence of external field. In Sec. IIVI 
we investigate the QD interaction with arbitrary state of 
quantum light in the weak driving field regime. Sect. IVlis 
devoted to the theoretical analysis of local field influence 



on the Rabi oscillations in the QD exposed to coherent 
states of light and Fock qubits. A discussion of the re- 
sults obtained is presented in Sect. |Vl] and concluding 
remarks are given in Sect. IVIII 

II. QUANTUM DOT - QUANTUM LIGHT 
INTERACTION: THEORETICAL MODEL 

A. Interaction Hamiltonian 

In this section we formulate the interaction Hamilto- 
nian for a QD exposed to quantized field accounting for 
the local-field correction. Later on we exploit the Hamil- 
tonian for the derivation of equations of motions describ- 
ing dynamical properties of this system. 

As aforementioned, the local field in QD differs from 
the incident one due to the d-d electrons-holes interac- 
tion. A general formalism accounting for the d-d interac- 
tions in atomic many-body systems exposed to photons 
has been developed in Refs. llslfT^ as applied to nonlin- 
ear optics of Bose-Einstein condensatesi^ii^. We extend 
this formalism to the case of the QD exciton driven by 
quantized light. 

Consider an isolated QD exposed to quantized electro- 
magnetic field. The electron-hole pairs in QD are as- 
sumed to be strongly confined; thus we neglect the static 
Coulomb interaction between electrons and holes. We de- 
compose the operator of the total electromagnetic field 
into two components. The first one, £v, represents a 
set of modes that do not contain real photons. The sec- 
ond component, £q, represents the set of modes emitted 
by the external source of light (real photons). Such a 
decomposition as well as the subsequent separate con- 
sideration of the field components is analogous to the 
Heisenberg-Langvein approach in the quantum theory of 
damping, see Ref. U. The total Hamiltonian of the sys- 
tem " QD-f electromagnetic field" is then represented as 

Ti = Ti.Q + Tiph + TYvac + Ti-w + ^iv , (1) 

where 7io,ph,vac are the Hamiltonians of the QD free 
charge carriers, the incident photons and the virtual vac- 
uum photons, respectively. The terms Tiio.iv describe the 
interaction of electron-hole pair with incident quantum 
field £q and with vacuum field £v, respectively. In the 
dipole approximation these Hamiltonians are given by 

^lo.iv ^-IJ^ C^^o,- + £o,.V)dh (2) 

where V is the QD volume and 'P{r,t) is the QD polar- 
ization operator. The Hamiltonian Tivac is as follows 

Hvac = ^ huJkvl^Vkx , (3) 
kX 

where vl.^ and Vkx are the creation and annihilation op- 
erators of vacuum photons, k is the mode index, indexes 
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A = 1,2 denote the field polarization. The operator of 
vacuum electromagnetic field is determined as 



(4) 



where f2 is the normalization volume and e^,;^ is the po- 
larization unit vector. 

As a first step in the development of our theory we 
eliminate from the consideration the exchange interac- 
tions and proceed to the direct interactions. For that 
purpose, we exclude the vacuum photon operators v^^x 
and v\,y^ expressing them (and corresponding Hamiltoni- 
ans Tivac and Tiiv) in terms of the polarization opera- 
tor 'p. Recalling the Heisenberg equation ifidvkx/dt = 



'H,vk\ 



the expression as follows can be obtained, 
dvkx 



dt 



-iLUkVkx + Fk\(t) 



where Fkx{t) = ^27rwfc/?il7 / -p{r' ,t)ekxe'^'^d^i 

V 

tion of Eq. (O is given by 



(5) 
Solu- 



Vkxit) = Vkx{-^)e~'^'* + / FfcA(r)e-^"-(*-^)dr (6) 



with the first term describing the free evolution of the 
reservoir modes (quantum noise) and the second one re- 
sponsible for the exchange interactions. Further we ne- 
glect the first term in ([5]) leaving the quantum noise be- 
yond the consideration. Inserting then this equation into 
(ID) and the resulting expression - into Hamiltonian 
after some algebra we arrive at 



2'KiLUk 



P„(r, t)ri3{r',T) 



-ooV V 



Y,^thfx(^'"^"~"'^('"""^"~'^dTSY<fv' + H.c. , (7) 



where indexes a, (3 mark Cartesian projections of vectors. 
The summation over repetitive indexes is assumed. Using 
the relationship^ 



E 



1 



fc2 \dxadxe, 



■ Sal 



1^ 



we proceed to the limit f2 oo in ([7]) . That corresponds 
to the replacement 



Then, utilizing the Markov property of the polarization 
operator, Va{'r,t) — Va{'r,0) the Hamiltonian ^ is 
reduced to 



-in 
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1 52 



vv 



. dxadxjs '"^ dt' 
X G^"'> {r - r' ,t')Va{r,t)Vp{r' ,t)dt'd^rd^r' , (8) 

where 

G(")(r,t) = 



2(27r)3 
1 

47r|r| 



zkr 



(e-^'^'=*-e*"'=*)d3k 



t 



(9) 



is the free-space Green function^ and 5{. . .) is the Dirac 
delta-function. Evaluation of Tivac in ^ is carried out 
analogously and gives Tivac = — '^iv/2. 

As the next step, we adopt the quasi-static approxima- 
tion, which utilizes the property of QD to be electrically 
small. The approximation implies the limit transition 
c cxD and the neglect the terms ~ /dt'"^ in Hamil- 
tonian ([8]). Then, the Hamiltonians Tivac and TLi^ are 
represented by the sum as follows 

A7i = Tivac + 'Hiv 

^ ^ ^'P(r)G(r-r')P(r')d^rd3r', (10) 



V V 



where 



G(r - r') = Vr «i Vr 



(11) 



is the free space Green tensor; Vr Vr is the operator 
dyadic acting on variables r. In the quasi-static approxi- 
mation we neglect the line broadening due to the dephas- 
ing and the spontaneous emission. The latter effect can 
be introduced in the model by retaining terms 0(l/c) in 
the the quasi-static approximation. 

In the preceding analysis we have suggested that the 
exchange by virtual photons of all modes occurs be- 
tween all allowed dipole transitions. That is, on that 
stage the problem was stated as a quantum-mechanical 
many-body problem. The analysis can be significantly 
simplified if we restrict ourselves to the two-level ap- 
proximation assuming the exciton transition frequency 
to be resonant with the acting field carrier frequency 
and utilize the self-consistent field model. The self- 
consistent field is introduced by means of the Hartree- 
Fock-Bogoliubov approximation's^ which implies the lin- 
earization of Hamiltonian (fTU|) by the substitution 

P(r)P(r') ^ V{r){V{v')) + {V{v))V[v') . (12) 
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The polarization operator of two-level system is given 



P(r) = |C(r)p(/xa+ + /x*<7_) 



(13) 



where a± are the Pauli pseudospin operators and C(r) is 
the wavefunction of the electron-hole pair. In the strong 
confinement regime this function is assumed to be turns 
out to be the same both in excited and ground states^i^. 

In the two-level approximation, the Hamiltonian of the 
carriers motion is represented as 



radius-vector of the QD geometrical center. The depo- 
larization tensor is given by 

K^-^IJ le(r)P |^(r')pG(r - r')d'rd^r' . (18) 

V V 

Noted that the resulting Hamiltonian (fT5|l coincides with 
that obtained in Ref. |35| in independent way. 



B. Equations of motions 



(14) 



where Sg^c and a| o/'^g.c ^-re the energy eigenvalues 
and creation/annihilation operators of the exciton; in- 
dices e and g correspond to the excitonic excited and 
ground states, respectively. The acting field operator 
is expressed by the relation ^ after the substitutions 



j(t) and 



5^ 



cjj{t); c]j{t)/cq are the cre- 



Vk\ - 

ation/annihilation operators of the incident (real) pho- 
tons (the polarization index A is included in the mode 
number q) . Formally, the relation ^ is fulfilled for oper- 
ators Cq{t) and cj(t) too, and the first term describes the 
evolution of real photons. However, since the exchange 
interaction is included into the vacuum field component, 
in the case of real photons the second term in relation 
([5]) disappears. Then, the Hamiltonian Tiph is given by 
the relation ([3]) after the substitution Vk\ — > Cq(— oo) and 



c*{—oo). For shortness, we denote Cq(— oo) 



and cj(— oo) = cj . 

Note that nonresonant transitions can be approx- 
imately accounted through a real-valued frequency- 
independent background dielectric function e^. Assum- 
ing eh to be equal to dielectric function of surrounding 
medium, we put further e/j = 1 without loss of gener- 
ality. Substitutions in final expressions c — > cj^fth and 
/X — + iij yflh for the speed of light and the electron-hole 
pair dipole moment, respectively, will restore the case 

As a next step we introduce the rotating wave 
approximation^, i.e., we neglect in ([8]) the terms that are 
responsible for the simultaneous creation/annihilation of 
exciton-exciton and exciton-photon pairs. Then, using 
expressions ([2]), ([3|) and ([8l)- (fT3)) . after some algebra we 
derive the effective two-particle Hamiltonian 



off 



where 



and 



^10 ^ f^^{9q^+Cq + 9*q^-cl) 



(15) 
(16) 



An 



An 
V 



fiiNf,){a^{a+)+a+{a.)). (17) 



where gq = ~ifj,eqy/2TTUJi^/hil exp(ikrc) is the coupling 
factor for photons and carriers in the QD and Tc is the 



Let \i^{t)) be a wavefunction of a QD interacting with 
quantum light. In the interaction representation the sys- 
tem is described by the Schrodinger equation 



ot 



(19) 



with \il){t)) = exp[i(7i;o + npx^)t/h]\'ip(t)) and Hint = 

exp[i(Ho + Hph)i/?i](^io + AH) exp[-(iHo + 'Hph)t/n]. 
We represent the wavefunction \ilj{t)) by the sum as fol- 
lows 

m))^ E [An.}(i)|e)+i?{„,}(t)|ff)] IK}), 

{«fc}>0 

(20) 

where A^n^y{t) and i3{„^}(t) are coefficients to be found, 
|{nfc}) denotes the multimode field state with n photons 
in k mode; |{Ofc)} is the wavefunction of the vacuum state 
of electromagnetic field; |e) and \g) are the wavefunctions 
of the QD ground and excited states, respectively. By 
inserting the relation (|20p into the Schrodinger equation 
(1191) . after some manipulations we arrive at the system 
of equations of motion 



dA 



{™i } 



dt 



dB 



dt 



+ 5]5,V^^^i?{,„,+5,je*('^»-"')* ,(21) 



with 



47r ~ 
Aw = ^M(iV/i) , 



(22) 



as the local-field induced depolarization shift^^'^^. Here 
^0 = (Eo — £g)/?i is exciton transition frequency. It can 
easily be shown that system (PT|) satisfies the conserva- 
tion law 



\B 



{n.k}\ 



= 0. 



(23) 
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The system ([21]) allows analyzing the interaction be- 
tween QD and electromagnetic field of an arbitrary spa- 
tial configuration and arbitrary polarization. Letting the 
coefficients gq{t) to be adiabatically slow-varying func- 
tions we can apply ((2T|) to QDs exposed to electromag- 
netic pulse. 

C. Single— mode approximation 

Among different physical situations described by Eqs. 
PT|) the single- mode excitation is of special interest. In- 
deed, such a case corresponds, for instance, to the hght 
-QD interaction in a microcavity with a particular mode 
resonant with the QD exciton. Owing to the high Q- 
factor, the strong light-QD coupling regime is feasible in 
microcavity providing numerous potential applications of 
such systemsii^. 

For the case of a spherical QD interacting with 
single-mode light, only the components |{?t-A;}) — 
|0i, O2, . . . riq . . . Ofe . . .) = \nq) with q as the number of 
interacting mode are accounted for in the wavcfunction 
(I20p . Then, omitting for shortness the mode number, the 
system (|22l) is reduced to 

+ 5V^rTTS„+ie^("»-")* , (24) 

■m 

Note that the conservation law ((23)) holds true for Eqs. 
([M)) with the substitution {rife} n. Equations (I^T]) 
and (IM)) govern the time evolution of the QD driven by 
quantum light. 

III. FREE MOTION 

The free motion regime implies neglect the QD- 
electromagnetic field interaction and thus imposes the 
condition g = on Eqs. (p4|) . Wavefunction of the nonin- 
teracting QD and electromagnetic field is factorized thus 
allowing analytical solution of ([M)) in the form of 

where C„ are arbitrary constants satisfying the normal- 
ization condition J2n — 1- that case the system 
(IM)) is reduced to the exactly integrable form 

^^^Ac.A\B\\ ^^ = Au;B\A\\ (25) 
and its solution is given by 

A{t) = aoe-^'^'^l''"!'* , B{t) = boe-'^'^^''"^'* . (26) 



Here ag and bo are arbitrary constants satisfying the con- 
dition |aoP -I- l^oP = 1- This solution describes a corre- 
lated motion of the electron-hole pair resulted from the 
local field-induced self-polarization of the QD. Thus, a 
quasi-particle with the wavefunction 

m)) = A{t)e-^'^^'\e) + Bit)e-^-^'\g) , (27) 

appears in the QD. It can easily be shown that the state 
(P7)) satisfies the energy and probability conservation 
laws. The inversion, which is defined as the difference be- 
tween the excited-state and the ground-state populations 
of the QD exciton, for the wavefunction (j27p remains 
constant in time: w = |^(i)P - \B{t)\'^ = |aoP - |6oP, 
whereas this state is generally non-stationary. The quasi- 
particle lifetime, which is not included in our model, can 
be estimated by Tgp 1/rsp, where Fsp is the QD- 
exciton spontaneous decay rate. For realistic QDs Aw ^ 
l/f^sp [111. Consequently, the state |i/'(0) can be treated 
as stationary within the range 1/Aw ^ t ^ 1/rsp- 

The macroscopic polarization of the QD is described 

by 

{V) = i^m) = l/xao&Se-^^"°-''^* + c.c. , (28) 

where the parameter 6' — wAoj plays the role of the self- 
induced detuning, which depends on the state occupied 
by the exciton and on the depolarization shift. Thus, as 
follows from ([28]), the local field-induced depolarization 
shift (Aw 0) dictates the non-isochronism of the polar- 
ization oscillations, i.e. the dependence of the oscillations 
frequency on its amplitude. This mechanism also influ- 
ences the Rabi oscillations in the system: the smaller 6' 
the larger Rabi oscillations amplitude; such a behavior 
was observed experimentally in Ref. uai 

Since the inversion w lies within the range — 1 < w < 1, 
the frequency ujp of polarization oscillations in (|28p may 
vary in the limits luq — Alo < LUp < loq + Alu. On the con- 
trary, when Alu — the polarization oscillates with the 
fixed frequency luq. At first glance, it seems that the dis- 
crete level is transformed into 2Auj band. However this 
is not the case. Indeed, the concept of the band structure 
corresponds to linear systems where any arbitrary state 
is a superposition of eigenmodes with different frequen- 
cies. As different from that, the electron-hole correlation 
arises from the nonlinear motion of the particles in a self- 
consistent field. Consequently, in the presence of light- 
QD interaction (g ^ 0), the exciton motion can not be 
described by a simple superposition of different partial 
solutions like (|26p . but has significantly more complicate 
behavior. In particular, in the strong coupling regime 
there exist two oscillatory regimes with drastically differ- 
ent characteristics separated by the bifurcation, see Sec. 

El 
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IV. THE WEAK-FIELD APPROXIMATION 

Consider a ground-state QD be exposed to an arbi- 
trary state of quantum light X^m, /^{mjli'^'})' where 
P{mi} 3.re arbitrary complex-valued coefficients satisfy- 
ing the condition X]{m,} l/^{ni;}P — 1- Then, the initial 
conditions for Eqs. (|2ip are given by 

A{„,}(0)=0, B{„,,}(0) (29) 

In the linear regime with respect to the electromag- 
netic field which realized when — > we can assume 
B{mi}{t) ~ P{m,} — const, i.e. the analysis is restricted 
to the time interval much less than the relaxation time of 
the given exciton state. Then, the system (1^^ is reduced 
to 



Now let us calculate the QD effective scattering cross- 
section defined by 



a(oo) = lim -^|a(t)p. 



(34) 



Substituting (j33|l into ((34|) , after some algebra we arrive 
at 



Cr((X)) = ^ ^ \C^q\'^d{uJ„ - UJq) , 



(35) 



dA 



dt 



- ^5,v/^^^/3{™,+.-„}e'^-°-')*.(30) 

For further analysis we rewrite (|30p in more convenient 
matrix notation: 



dSL{t) 

dt 



= -iAujpa.{t) - f (t) 



(31) 



where Ciq = [p — I)fq and = pfq, loi = loq. The 
equation obtained demonstrates a fine structure of the 
absorption in a QD exposed to quantum light with arbi- 
trary statistics. Instead of a single line with the exciton 
transition frequency wi = wq, a doublet appears with one 
component shifted to the blue: u)2 = ujq + Aw. The in- 
tensities of components are completely determined by the 
quantum light statistics. The same peculiarity has been 
revealed in Ref. [s^ in a particular case of a QD exposed 
to the single-mode Fock qubit. Obviously, this result di- 
rectly follows from psp if we retain in this equation only 
terms corresponding the mode considered. The single- 
mode Fock qubit is a superposition of two arbitrary Fock 
states with fixed number of photons and is described by 
the wavcfunction (the mode index q is omitted) 



where a(i), i{t) = f^e^^""""")* and /3 are the colum- 
nar matrices (vectors) and a(i) = (A{mj). Vectors 
fq and (3 are defined analogously through the elements 

/j'"'^ = igq^/rug + l/3|,„,+5^j and /3{™j, respectively; 
6qi is the Kroneker symbol. The quantity p — f3f3^ is the 
density matrix of quantum light interacting with the QD. 
This matrix corresponds to the pure state of electromag- 
netic field and satisfies the conditions p" — p {n > 1) 
and I (where I is the unit matrix). Integration of 

([3T|) imposed by initial conditions ([29]) allows us to find 
the analytical solution 



^f3N\N) + l3N+i\N + l). 



(36) 



of 



a(i) = - y e-'P'^'^'''-*'k{t')dt' 



Using the truncated Taylor expansion^ 



(32) 



-ipAuj^t — t') 



after some standard manipulations we obtain 



a(i) 



I)f„ 



1 



piqC 



^ LOii — LOq — iQ 
_ gi(wo-i^<j+Aii;)ty 
UJ(i — UJq + Aw — iO 



(33) 



This expression describes the interaction between the QD 
and an arbitrary state of quantum field in the weak-field 
limit. 



Consequently, at > 1 the nonzero components 
the vector f, are /(^-i) = igVNpN and /(^) = 
ig^/N + IPn+i ■ It can easily be found that in that case 
Eq. ((35|) is reduced to Eq. (61) from Ref. 35. In the case 

= the only nonzero component J*^*^' = wPi survives 
in (|35p thus reducing this equation to Eq. (67) from Ref. 

For the single Fock state |-0) — \N) and we obtain 
C2 — 0, i.e., in agreement with Ref. [H, only shifted 
spectral line is presented in the effective cross-section. 

For the QD ex posed to coh erent states the recurrent 
formula /3n+i = a/ {n)/{n + l)/3„ can easily be obtained, 
where (n) stands for the photon number mean value. Us- 
ing this formula the relation pi = f can be obtained, 
which gives Ci = 0. Thus, for this case only shifted line 
in the effective cross section is manifested. Analogous re- 
sult has been reported in [3^ for a QD driven by classical 
electromagnetic field. Such a coincidence can be treated 
as the manifestation of well known concept: among a va- 
riety of quantum states of light the coherent states are 
most close to classical electromagnetic field. 



V. THE QUANTUM LIGHT-QD STRONG 
COUPLING REGIME. RABI OSCILLATIONS 

In the case when the coupling constant is compara- 
ble in value with the exciton decay rate, the linearization 
with respect to electromagnetic field utilized in the previ- 
ous section is no longer admissible. The strong coupling 
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regime can be realized by combining the QD with a high- 
Q microcavity and one of its manifestation is the Rabi 
oscillations of the levels population in a two-level sys- 
tem exposed to a strong electromagnetic wave. In this 
section we utilize the system of equations (|24l) for inves- 
tigation of the interaction between QD and single-mode 
quantum light. The incident field statistics is exempli- 
fied by coherent states and Fock qubits. In the standard 
JC model the Rabi oscillations picture is characterized 
by the following parameters-i: (i) the coupling constant 
g] (ii) the frequency detuning 5 ~ lo^) — lo] (iii) the initial 
distribution of photonic states in quantum light and (iv) 
the average photon number {n). Accounting for the local 
field correction supplements this set with the new pa- 
rameter, the depolarization shift Aw. For convenience, 
we introduce the depolarization shift by means of the 
parameter-- ^ = r2(„^ /Acj which compares the shift with 

the Rabi frequency — 2|g|y^n(0)). To understand 
the dynamical property of the Rabi effect we have inves- 
tigated the following four physical characteristics: the 
QD-exciton inversion 



;(t) = ^(|A„(t)p-|i?„(t)r), 



(37) 
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the time evolution of the distribution of photonic states 

p(n,t) = |A„(t)p + |B„(t)|2, (38) 

and the normally ordered normalized time-zero second 
order correlation function of the driving field: 



= J2nin-l)pin,t)/[Y^np{n,t)]\ (39) 



The temporal evolution of the inversion can be de- 
tected in pump-probe experiments^ii^, while the pho- 
tonic state distribution psp is measurable in quantum 
nondemolition experiments with atoms^. 



A. Coherent states excitation 

Let a ground-state QD be exposed to the elementary 
coherent state of light X^J^o where-^ F{n) = 

exp[-(n(0))/2](n(0))"/Vv^- Then, initial conditions 
for Eqs. ([M]) are given by 



B„(0) = Fin) 



A.(0) = 0. 



(40) 



Figures [T] and [5] show calculations of the inversion in 
a lossless QD as a function of the dimensionless time 
T = \g\t at the exact synchronism {S = 0) for two different 
initial photon mean numbers {n(0)) and for several val- 
ues of parameter ^. Our calculations demonstrate the ap- 
pearance of two completely different oscillatory regimes 



FIG. 1: Rabi oscillations of the inversion for a QD exposed 
to the coherent state of light with (n(0)) = 9 and ^ = 0.2 (a), 
^ = 0.49 (b), C = 0.53 (c), C = 1.2 (d), C = 3.5 (e), C = 18.0 
(f)- 



in the Rabi oscillations. The first one manifests itself 
at ^ < 0.5 and is characterized by periodic oscillations 
of the inversion within the range — 1 < w{t) < (see 
Figs. [1^, [H)). Thus, in this regime the inverted pop- 
ulation is unreachable. On the contrary, in the second 
regime, at ^ > 0.5, the inversion oscillates in the range 
-1 < w{t) < 1 (Figs. Ht-E and Fig. [J). These two 
regimes of the Rabi effect are separated by the bifurca- 
tion which occurs at ^ = 0.5 for both types of incident 
coherent states (compare Figs. [T]d and[TJ;). 

In the limit ^ — *■ oo (Aw — > 0) the contribution of 
terms 0{Auj) in Eqs. is small. The neglect these 

terms corresponds to the elimination of the local field 
effect. In this case the system (1^^ is reduced to that 
follows from the standard JC model^, allowing thus the 
analytical solution: 



wit) 



n=0 



cos(f2„t) 



(41) 



where f2„ = ^/W+A\g^^^in + Y) . The fundamental effect 
predicted by this solution is the collapse-revivals phe- 
nomenon in the time evolution of the inversion^ . We have 
found that at ^ > 18 the numerical simulation by Eqs. 

leads to the same result as analytics (|H|), see Fig. 
[Tf. In the case ^ the amplitude of Rabi oscillations 
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The collapse-revivals effect in the time evolution of the 
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FIG. 2: Rabi oscillations of the inversion for a QD exposed 
to the coherent state of light with (n(0)) — 36 and ^ — 0.53 
(a), ^ = 1.2 (b), C = 18.0 (c). 
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FIG. 3; The second order time-zero correlation function 
g^^\t) of QD exposed to the coherent state of hght with 
{n(0))=9 and ^ = 18.0 (1), ^ = 3.5 (2), ^ = 1.2 (3) 

For a single QD imposed to classical light the appear- 
ance of two oscillatory regimes in Rabi oscillations sep- 
arated by the bifurcation at £ = 0.5 has been predicted 
in Ref. [Sy. Accordingly to [30|, the region ^ > 0.5 corre- 
sponds to periodic anharmonic oscillations of the inver- 
sion. As different from that, in a QD exposed to quan- 
tum light the collapse-revivals phenomenon takes place 
in this region, see Figs. [TJ;-f and[2l As Figs. [Tf and[2j; 
demonstrate, the collapse and the revivals in the vicinity 
of the bifurcation are deformed and turn out to be drasti- 
cally different from those predicted by the solution (|¥T|) . 



FIG. 4: Photonic state distribution in a QD exposed to co- 
herent light with (n(0)) = 9 for ^ = 1.2 and r = (a), r = 8.5 
(b), r = 14.5 (c), r = 31 (d), r = 40 (e). 

inversion disappears completely in the range ^ < 0.5 (see 
Figs. H^-b), where the Rabi effect picture turns out to 
be identical to the case of QD excited by classical light 

Let us estimate material parameters which provide ob- 
servability of the effects predicted. For a spherical In- 
GaAs QD with 6 nm radius the dipole moment can be 
estimated^^ as /x « 12 Debye. For this QD we obtain 
hAui « 0.1 meV. Then, for the range of ^ presented in 



Fig. m from f = 0.2 to ,f = 18, we obtain hn 



and hi} 



(n) 



0.02 



2 meV, respectively. These values are of 



the same order as the excitonic Rabi splitting measured 
in recent single QD spectroscopy exper iments, see Refs. 
[miii. On the other hand, Refs. SliillOIEl! report the 
QD exciton linewidth Fhom of the order of 1 /ieV below 
the temperature 10 K and laying in the range 4 to 10 
/ieV at T = 20 K. Thus, the precondition to observe the 
strong couphng regime, fi^{n) ^ rhom, is fulfilled for the 
given range of ^. 

An important feature of the Rabi effect in quantum 
light is the variation of the light statistics during the 
interaction with a quantum oscillator (QD). We shall 
characterize the variation by the second-order time-zero 
correlation function and the photonic state distribution 
defined by Eqs. pQ]) and ((38)) . respectively. These char- 
acteristics at different ^ are depicted in Figs. [3] and U) 
One can see that at large ^ > 18 the function g^'^\t) os- 
cillates around unit in agreement with the standard JC 
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modelA, see curve 1 in Fig. [3l The situation is changed in 
the vicinity of the bifurcation as it is presented by curves 
2 and 3 in that figure. As one can see, at ^ = f.2 the 
hght statistics becomes super-Poisonnian. The correla- 
tion function g^'^\t) demonstrates the increase in time 
imposed to small-amphtude oscillations. These oscilla- 
tions correspond to regions of revivals in the time evo- 
lution of the inversion (compare curve 3 in Fig. [3] and 
Fig. [T]i). It should be noted that below the bifurcation 
threshold, at ^ < 0.5, g'^'^\t) oscillates in the vicinity of 
unity too. 

Fig. m presents the photonic state distribution p{n, t) 
for different time points at the given ^ = 1.2. The 
figure illustrates consecutive transformation of the ini- 
tially Poissonian distribution (Fig. [4^) into the super- 
Poisonnian one in the course of time, see Figs. [31;-d. 
The transformation corresponds to the increase in g*^^^ (t) 
illustrated by curve 3 in Fig. [3l Let us stress that the 
standard JC model predicts the photon statistics remain- 
ing Poisson as it takes place in our case only at large ^ 
(curve 1 in Fig. [3]). Our calculation also show that in 
the region ^ < 0.5 the photon distribution p{n,t) does 
not depend on time and remains Poissonian. The invari- 
ability of the coherent light statistics in the limit ^ ^ 
corresponds to the absence of the component h' = 1 in 
Eq. dSSI). 
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FIG. 5: Rabi oscillations of the inversion for the QD exposed 
to the Fock qubit for ^ = 1.2 (a), ^ = 8.0 (b), ^ = 55.0 (c). 



B. QD interaction with Fock qubits 

Let a ground-state QD interacts with electromagnetic 
field given by the Fock qubit (f36|) . Then the initial con- 
ditions for Eqs. (f24|) are given by 

A„(0) = 0, Bn{0)-^SnMl3N +SnM + l/3N+l, (42) 

Further we restrict ourselves to t he Fock qubit with N — 
6 assuming (3n = Pn+i = \/l/2. 
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FIG. 6: Second order time-zero correlation function g^^\t) of 
the QD exposed to the Fock qubit with = 6 and ^ = 55.0 
(i),^ = 8.0 (2),C = 1.2 (3). 



Calculations of time evolution of the inversion for dif- 
ferent values of ^ are shown in Fig. [51 At large ^ 
(Fig. [5j:) the local field effect is eliminated: In agree- 
ment with the solution (|4ip the inversion exhibits har- 
monic modulation of the oscillation amplitude within the 
range we [—1,1]- The oscillation frequency is equal to 
(JIat -|- flN+i)/'2 while the frequency of the modulation is 
given by {Hn — ^n+i)/2. With the parameter ^ decrease 
the modulation becomes non-harmonic. 

The next two figures, [S] and [71 illustrate the change of 
the light statistics due to interaction with QDs. The time 
evolution of the correlation function g^^-* (t) is depicted in 
Fiq. [HI for different values of Curve 1 demonstrates 
that, same as in Fig. [31 and in agreement with the stan- 
dard JC model, at large C > 55) the correlation func- 
tion oscillates in the vicinity of its initial value. Curve 
3 shows transformation of the sub-Poissonian statistics 
{g^'^H^) < 1) info the super-Poissonian one (g'^^^t) > 1) 
with a pronounced maximum followed by irregular oscil- 
lations. Analogous but essentially smoothed behavior is 
observed at larger ^ (curve 2). 

Figure \7\ presents calculations of the photonic state 
distribution p{n,t) for ^ = 1.2 (curve 3 in Fig. [5]) in 
different points of time. The standard JC model when 
describes the interaction of Fock qubit with two-level 
system, predicts variation of those probabilities p{n, t) 
([55)1 that correspond to Fock states with n = TV, ± 1 
numbers of photons. As different from that, the incorpo- 
ration of the local-field effect leads to the appearance in 
the distribution of states \n) with photon numbers both 
smaller and larger than present in the initial Fock qubit 
(Fig[7^). Probabilities of these states are redistributed 
(Figs. Wp^c) with time, and light statistics became irreg- 
ular which signifies the transformation of the photonic 
state distribution (Figs. [71i,e). In turn this affects the 
Rabi oscillations picture and the second order correlation 
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FIG. 7; Photonic state distribution of the QD exposed to the 
Fock qubit for ^ = 1.2 and r = (a), r = 5.19 (b), r = 14.5 
(c), T = 29.6 (d), T = 40.1 (e). 



function g'''^\t) (compare Figs. [5^ and[6l curve 3). 



C. Vacuum Rabi oscillations 

The vacuum Rabi oscillations characterize interaction 
of an excited-state QD with electromagnetic vacuum. 
The initial conditions for Eqs. l[24l) in that case are given 
by Bn{Q) = 0, ^n(O) = 5n,o- Numerical solution of this 
system leads to the time-harmonic oscillations of the in- 
version. This agrees with the analytical solution of the 
system at Aw — given by the standard sinusoidal law^: 



4|aP 

cos(floi) . 

"o 



(43) 



The result can easily be understood from that the vac- 
uum states, like a single Fock state, have zero observ- 
able electric field. Therefore, such states do not induce 
in QDs observable polarization and, consequently, fre- 
quency shift. For the same reason, zero frequency shift is 
inherent to single-photon states, as it has been revealed 
under some simplifying assumptions in Ref. [ssl . 



VI. DISCUSSION 

The standard JC model describing the interaction be- 
tween single-mode quantum electromagnetic field and 
two-level system predicts the collapse-revivals picture of 
the time evolution of the level population. The basic 
physical result of the analysis presented in our paper is 
a significant modification of the Rabi effect due to the 
local-field induced depolarization of a QD imposed to 
quantum light. Two oscillatory regimes with drastically 
different characteristics arise. In the first regime the time 
modulation of the population (the collapse and revivals) 
is suppressed and the QD population inversion is found 
negative. This indicates that trajectories of charge car- 
riers confined by the QD occupy a finite volume in the 
phase space. In the second oscillatory regime the re- 
vivals appear, however they are found deformed and sig- 
nificantly different from that predicted by the standard 
JC model. The trajectories of charge carriers occupy the 
entire phase space. Both regimes of oscillations demon- 
strate the non-isochronous dependence on the coherent 
field strength. 

Two regimes of the Rabi oscillations indicate the ap- 
pearance of two types of motion of the QD exciton. The 
first one is a superposition of time-harmonic oscillations 
with the Rabi frequencies il„ = g^N + 1, while the sec- 
ond type is presented by the frequency band \ujp — o^ol < 
2Aw . The resulting exciton motion is thus determined 
by a nonlinear superposition of these two types of motion. 

The first mechanism of the motion is conventional for 
the Rabi effect; physically it originates from dressing of 
the QD exciton by the incident field photons. This type 
of motion dominates at large strengths of the incident 
field, when ^[n) ^ Aw. With the light-QD coupling 
constant (and correspondingly r2(„)) decrease, the role 
of the second type of motion grows in importance. The 
reason for this is the electron-hole correlations resulting 
from the exchange interaction. It can be interpreted as 
the QD exciton dressing by virtual photons. This regime 
becomes dominating in the comparatively weak fields, 
when Vli^n) ^ ^i^- Thus, the reduction of the threshold of 
the acting field strength needed for the Rabi oscillations 
appearance, recently observed experimental!}^, can be 
attributed as a local-field effect. The experimen1>2^ has 
also elucidated the non-isochronism of excitonic Rabi os- 
cillations that can be treated as the local-field effect as 
well; see the relation (l28l) and numerical calculations re- 
ported in Ref. [36]. 

It should be noted that two oscillatory regimes in the 
Rabi effect may appear in other quantum systems where 
additional interaction mechanisms exist. As an exam- 
ple, consider a two-component Bose-Einstein conden- 
sate with radio- frequency coupling of two separate hyper- 
fine states^. Temporal evolution of this system is gov- 
erned by the coupled Gross-Pitaevskii equations, which 
are similar to those derived in Sec. IIIBI The equations 
combine both linear and nonlinear couplings. The linear 
coupling constant characterizes the interaction between 
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the system and the electric field, while the nonhnear one 
accounts for the interaction between intra- and inter- 
species of the condensate^. Dependently on ratio of the 
couphng parameters, the Rabi oscillations between the 
condensate components may exhibit both well-ordered 
and chaotic behavior^^, similar to that depicted in Fig. 
[T^,b. The formation of the Bardeen-Cooper-Schrieffer 
state in the fermionic alkali gases cooled bellow degener- 
acy [l^l can serve as another example. In that system, 
the time modulation of coupling constant leads to the 
Rabi oscillations of the energy gap^ with two oscillatory 
regimes. The trajectories of individual Cooper pairs oc- 
cupy a finite volume in the phase space in the first regime 
and the entire phase space in the second one. 

Now, let us turn to the weak-field case. The the de- 
polarization induced local field is predicted to entail in a 
QD exposed to an arbitrary photonic state a fine struc- 
ture of the effective scattering cross-section. Instead of a 
single line of the frequency loq , a duplet is appeared with 
one component shifted by a value Aw. The shifted com- 
ponent is due to electron-hole correlations, see Eq. 
The correlations change the QD state and, consequently, 
provide the inelastic channel of the light scattering. The 
elastic scattering channel is formed by light states induc- 
ing zero observable polarization and, consequently, zero 
frequency shift, such as Fock states, vacuum states, etc. 
Now we take into account the relation (P) = 47rQ!(Eo), 
which couples observable polarization and mean value 
of the incident field; the quantity Eq is defined through 
the relation 5o(t) — Eo(r, ti;)(ia; -I- H.c. The scalar 
coefficient a is the QD polarizability of a spherical QD. 
Therefore, we conclude that the elastic scattering channel 
is formed by incident field with zero mean value (incoher- 
ent component of the electromagnetic field). Correspond- 
ingly, the coherent field component is scattered through 
inelastic channel. As follows from the solution ((35| . the 
elastic channel is not manifested for pure coherent light 
(the Glauber state). 

Let us discuss now the local-field induced alteration of 
the quantum light statistics. As an example we consider 
the Fock qubit ([M)) as the incident field state. For the 
case Alli = the photonic state distribution is given by 



p{n,t) = Sn,N^l\AN-lit)\^ +Sn^N{\AN{t)\'' 



(44) 



where An{t) and Bn{t) are the exact solutions of Eqs. 
(1^ at Au) — 0, see, e.g. Ref. HI It is seen that the Fock 
states with photon numbers n = N, N ± 1 are only pre- 
sented in the distribution. The probability amplitudes 
of these states oscillate with the corresponding Rabi fre- 
quencies, ilii=N.N±i- In addition to this set, extra Fock 
states with both smaller and larger photon numbers are 
appeared in the photonic state distribution p{n, t) as the 
local-field effect. See Fig. [T] Therefore, a bigger number 
of Fock states than presented in the initial Fock qubit, 
broaden the frequency spectrum of Rabi oscillations, pro- 
viding thus the chaotic time evolution of the inversion. 



It should be noted that the variation in the quantum 
light statistics occurs even in the weak-field limit 5 — > 0. 
To illustrate that, consider the observable polarization of 
the QD exciton defined in the time domain as 



1 
V 



E 



A, 



+ C.C . 



(45) 



Using Eq. ([33)1 we couple the polarization in the fre- 
quency domain with the complex-valued amplitude of 
the mean incident field: 



ujQ ~ uj + Auj — iO 



(En 



(46) 



Then, after some simple manipulations we express the 
quantum fluctuations of the QD polarization by 



P (P) 



ujo — bj — iQ 



Eo — (Eo' 



(47) 



From Eqs. (gej) and follows that the QD electro- 
magnetic responses to the mean electric field and to its 
quantum fluctuation are different: the response resonant 
frequency is shifted by the value Auo in the former case 
[relation (|46|) ] and remains unshifted in the later one, 
as given by relation ([17|) . This indicates that the effec- 
tive polarizability of a QD is an operator in the space of 
quantum states of light. It should be pointed out that 
this property is responsible for the alteration of the pho- 
tonic state distribution in the weak-field regime and is 
entirely a local-field effect. Note that the notions "strong 
(weak) coupling regime" and " strong (weak) field regime" 
are not identical as applied to QDs. To illustrate this 
statement, we express from Eqs. and (HT)) the po- 

larization operator in the weak-field limit: 



ujQ — uj — i{) \ 



AujCEn 



uja — Lo + Auj — iO 



(48) 



This equation is linear in the incident field but includes a 
term quadratic in the oscillator strength, 0(|/x|'*). Non- 
linearity of that type violates the weak coupling regime. 
The splitting of the QD-exciton spectral line dictated by 
Eq. (|47)l is a manifestation of the strong light-matter 
coupling in the weak incident-field regime. Thus, the 
light-QD interaction is characterized by two coupling pa- 
rameters, the standard Rabi-frequency and a new one, 
the depolarization shift Au. 



VII. CONCLUSIONS 

In the paper we have developed a theory of the elec- 
tromagnetic response of a single QD exposed to quantum 
light, corrected to the local-field effects. The theory ex- 
ploits the two-level model of QD with both linear and 
nonlinear coupling of excited and ground states. The 
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nonlinear coupling is provided by the local field influ- 
ence. Based on the microscopic Hamiltonian accounting 
for the electron-hole exchange interaction, an effective 
two-body Hamiltonian has been derived and expressed 
in terms of the incident electric field, with a separate 
term responsible for the local-field impact. The quan- 
tum equations of motion have been formulated and solved 
with this Hamiltonian for different types of the QD exci- 
tation, such as Fock qubit, coherent state, vacuum state 
and arbitrary state of quantum light. 

For a QD exposed to coherent light we predict two os- 
cillatory regimes in the Rabi oscillations separated by the 
bifurcation. In the first regime, the standard collapse- 
revivals phenomenon do not reveal itself and the QD 
inversion is found negative. In the second regime, the 
collapse-revivals picture is found to be strongly dis- 
torted as compared with that predicted by the stan- 
dard JC model. The model developed can easily be ex- 
tended to systems of other physical nature exposed to 
a strong electromagnetic excitation. In particular, we 
expect manifestation of the local-field effects in Bose- 
Einstein condensates^ and fermionic alkai gases cooled 
below the degeneracj^'^. 

We have also demonstrated that the local-field correc- 
tion alters the light statistics even in the weak-field limit. 
This is because the local fields give rise to the inelastic 
scattering channel for the coherent light component. As 
a result, coherent and incoherent light components in- 
teract with QD on different frequencies separated by the 
depolarization shift Aw. In other words, the local fields 
eliminate the frequency degeneracy between these com- 
ponents of the incident light. 

Note that our model does not account for the 
dephasing. Accordingly to recent experimental 
measurements^!^ and theoretical estimate o^^i^^ , the 
electron-phonon interaction is the dominant mechanism 
of the dephasing in QDs. Thus, the further development 
of the theory presented requires this dephasing mecha- 
nism incorporation. A next step is generalization of our 
model to multi-level systems. Among them, the systems 
with dark excitons interacting with weak probe pulse in 
the self-consistent transparency regimeil are of special 
interest. 

In the paper we have considered an isolated QD. The 
generalization of the theory developed to the case of QD 
ensembles (excitonic composites), such as self-organized 
lattices of ordered QD molecules^ and ID-ordered (In, 
Ga)As QD arrays^, is of special interest. One can expect 
that dipole-dipole interactions between QDs will mani- 
fest itself in a periodic transfer of excited state between 
QDs resulting thus in the collective Rabi oscillations — 



Rabi waves. 
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APPENDIX A: DEPOLARIZATION SHIFT FOR 
THE SPHERICAL QD 

For a spherical QD, the formulae is reduced to 

Au; = ^lMpTr(iV). (Al) 

Using the expression p?]) we obtain 

Tr(^) = /|C(r')p/lC(r)pV^(^)dW 

V V 

= V J y"|C(r')|'C(r)l''5(r-r')d^rd3r' 

V V 

= vj \C{r)\'d^r. (A2) 

V 

Consider the Is state. The normalized wavefunction in 
this case is given by^ 

where J^{. . .) is the Bessel function, p is the spherical 
coordinate and R is the QD radius. Inserting (jA3p into 
(IA2p and integrating the resulting expression, we derive 
the correction to the depolarization shift (|Aip 

„ ,~ ^ 47r [ sin x , . , 

Tr(iV = — / — —dx « 2.81 . A4 
3 J x^ 
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